Electron Momentum Distribution Function in the t-t'-J Model 

A. Ramsak^'^, I. Sega^, and P. Prclovsek^'^, 
J. Stefan Institute, 1000 Ljubljana, Slovenia 
^ Faculty of Mathematics and Physics, University of Ljubljana, 1000 Ljubljana, Slovenia 

(February 1, 2008) 






o 

(N 



We study the electron momentum distribution function (EMDF) for the two-dimensional t-t'-J 
model doped with one hole on finite clusters by the method of twisted boundary conditions. The 
results quantitatively agree with our analytical results for a single hole in the antiferromagnetic 
background, based on the self-consistent Born approximation (SCBA). Moreover, within the SCBA 
an anomalous momentum dependence of EMDF is found, pointing to an emerging large Fermi sur- 
face. The analysis shows that the presence of next-nearest-neighbor (NNN) hopping terms changes 
EMDF only quantitatively if the ground state (GS) momentum is at (f , f ) and qualitatively if the 
GS momentum is shifted to (vr, 0). 
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I. INTRODUCTION 

The existence and the character of the Fermi surface 
(FS) of high Tc superconductors, in particular in their 
underdoped regime, is one of the open questions ||l[ of 
solid state physics. The key quantity for resolving this 
problem is the electron momentum distribution function 
"k = (*kol E<T 4,<TCk,<T|*ko)- Here we study the EMDF 
for I'I'ko) which represents a weakly doped antiferromag- 
nct (AFM), i.e., it is the GS wave function of a planar 
AFM with one hole and the GS wave vector kp. The wave 
function is determined within the framework of the stan- 
dard t-J model with nearest-neighbor hopping tjji = t 
and the AFM exchange is fixed to J = 0.3t. In order 
to come closer to the realistic situation in cuprates the 
model is extended with the NNN hopping tjji = t' for jj' 
representing next-nearest-neighbors, 



H=- Y. t 

<jj'>c 



n' y'^j.a^j' ,<r 



H.c. 



<u> 



+J ^ [StS^ + ^{StSJ + Si 8+)] . 



(1) 



Numerically rik can be determined by exact diagonaliza- 
tion (ED) of the model Eq. (|l|) in small clusters, where 
only a restricted number of momenta k is allowed. The 
GS wave vector due to finite size effects varies with 
N. Therefore we present here results obtained with 
the method of twisted boundary conditions B, where 
tjji — > tjji expiOjji. Since rik = nk(ko,0) depends both 
on kg and 9 it follows from Peierls construction that 
nk(ko,0) = n-k+ko(Oiko) for 9 = kg. This allows us to 
study TT-k for arbitrary k and kg. Furthermore, the finite 
size effects of the results are suppressed if we fix kg for all 
clusters here studied to the symmetry point ko = (f-, §)■ 



II. RESULTS 

In Fig. 1(a) we present ED results for clusters with 
different N and 7 = 1. The EMDF obeys the sum 
rule ^j^ rik = A^ — 1 and, for the allowed momenta, the 
constraint N{n\t, — 1) < 1. We show here the quan- 
tity N(nu — 1), which for different N scales towards 
the same curve. Results are presented for momenta 
k e [(— TT, —n) -^ (tt, tt)] and should be averaged over 
all four possible ground state momenta when discussed, 
e.g., in connection with ARPES data. 

We further compare the results with the self-consistent 
Born approximation M, where we decouple fermion op- 
erators into hole and magnon operators |g|,g[, and us- 
ing the SCBA wave function |H,0 we evaluate the re- 
quired matrix elements in nu- In Fig. 1(a) the result for 
N{nu — 1) is presented. In the SCBA there appear in n^ 
at k ~ ko,ko -I- (7r,7r) also delta- function contributions 
proportional to the quasiparticle pole residue, Zko, but 
not shown in the figures. The comparison of the SCBA 
with the ED results shows a quantitative agreement, how- 
ever, a surprising observation is that n\^ in the SCBA 
exhibits a discontinuity for momenta k ~ kg, ko -I- (tt, tt) 
m . We interpret this result as an indication of an emerg- 
ing large Fermi surface at k ~ ikg indicatingt the co- 
existence of two apparently contradicting FS scenarios 
in EMDF of a single hole in an AFM. On one hand, 
the (5-function contributions seem to indicate that at fi- 
nite doping a delta-function might develop into small FS 
(hole pockets). On the other hand the discontinuities at 
k = ko, ko + (tt, tt) are more consistent with infinitesi- 
mally short arcs (two points) of an emerging large FS. 

In Fig. 1(b) and Fig. 1(c) we also present the results 
for the t-t'-J model, (i) The effect of negative NNN t' = 
— i/4 is relatively weak: the GS momentum remains at 
ko — (f J f )■ Otherwise results are similar to the t' — 
case, (ii) Positive NNN hopping matrix elements t' — 
t/4: the GS is now twofold degenerate, with the momenta 



at corners of the AFM zone, e.g., ko = (7r,0). All results 
agree with the SCBA method. 
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FIG. 1. N{n\^ — 1) obtained from ED for various systems 
N = 16, 18,20 and Af = 32 from Ref. [3]. For the SCBA 
Af < 64 X 64, 7 = 0.999 and note that delta-function contri- 
butions at k = ±ko are not shown, (a) t-J model, t' — 0, (b) 
t-t'-J model with t' — ±t/4 for k G [(— tt, — tt) -^ (tt, n)], and 
(c) for k G [(— vr, tt) -^ (tt, —it)]. 



